Absorption, Photoluminescence and Resonant Rayleigh Scattering Probes of 

Condensed Microcavity Polaritons 
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We investigate and compare different optical probes of a condensed state of microcavity polari- 
tons in expected experimental conditions of non-resonant pumping. We show that the energy- and 
momentum-resolved resonant Rayleigh signal provide a distinctive probe of condensation as com- 
pared to, e.g., photoluminescence emission. In particular, the presence of a collective sound mode 
both above and below the chemical potential can be observed, as well as features directly related 
to the density of states of particle-hole like excitations. Both resonant Rayleigh response and the 
absorption and photoluminescence, are affected by the presence of quantum well disorder, which 
introduces a distribution of oscillator strengths between quantum well excitons at a given energy and 
cavity photons at a given momentum. As we show, this distribution makes it important that in the 
condensed regime, scattering by disorder is taken into account to all orders. We show that, in the 
low density linear limit, this approach correctly describes inhomogeneous broadening of polaritons. 
In addition, in this limit, we extract a linear blue-shift of the lower polariton versus density, with a 
coefficient determined by temperature and by a characteristic disorder length. 

PACS numbers: 71.36.+C, 78.35.+C, 78.40.Pg, 71.45.-d 



I. INTRODUCTION 

Since the prediction of Kcldysh and Kopaevji there 
has been a long and intense pursuit to realise a con- 
densed phase in solid state excitonic, and related, sys- 
tems. In particular, polaritons in semiconductor micro- 
cavities, the coupled eigenstates of an exciton with a 
cavity photon^^ represent ideal candidates for observing 
condensation phenomena. The very light mass of these 
composite bosonic particles promises relatively high tran- 
sition temperatures. In the past decade, improvements in 
the growth technology of semiconductor heterostructures 
have made the study of high-quality strongly-coupled pla- 
nar microcavities almost routine for III-V and II- VI semi- 
conductors. The high degree of external control of these 
systems, and the possibility of direct detection for them 
has opened the route towards a new generation of fast op- 
tical matter-wave lasers and amplifiers.— More re- 
cently, concerted experimental efforts have been devoted 
to the realisation of a Bose-Einstcin condensate of micro- 
cavity polaritons^i 10 i u i 12 ' 13 i 14 

On the experimental side, a challenge to the reali- 
sation of a condensed polariton phase might be repre- 
sented by the finite quality of the cavity mirrors and 
the resultant short polariton lifetime, of the order of pi- 
coseconds. In addition, due to the 'bottleneck effect',— 
the relaxation of polaritons to the zero momentum state 
can be delayed, hindering the creation of a thermal 
population in the lowest energy state. It has however 
been recently show n 13 ' 14 ! 16 that thermalization processes 
due to particle-particle scattering can be dramatically 
magnified by increasing the value of the (non-resonant) 
pump power, and by positively detuning the cavity en- 
ergy above the excitonic energy. Under these conditions, 



the progress towards a zero momentum quasi-equilibrium 
condensate has been significant j 4 ^! 9 - ' 10 ' 11 ' 12 ' 13 including a 
non-linear threshold behaviour in the emission intensity 
at zero momentum,— the investigation of the sec- 
ond order coherence function^ 9 - a characteristic change in 
the momentum space distribution above threshold ) 10 ' 13 
and evidence that the equilibration time is shorter than 
the polariton lifetime has been seen.— Finally, very re- 
cently, a clear demonstration of condensation of cavity 
polaritons has been demonstrated in CdTe.— Kasprzak 
and collaborators have shown that condensation of equi- 
librated polaritons can be achieved for effective tempera- 
tures around 20K and evidence for condensation has been 
seen in the occupation function, the first order coherence 
(both in time and in space) and in the spontaneous ap- 
pearance of linear polarisation of the condensate emis- 
sion. 

Alongside the experimental effort, a significant the- 
oretical effort has been invested in analysing proper- 
ties and predicting signatures of polariton condensa- 
tion, 1 ^ 18 ! 19 . 20 ! 21 ^ 2 ! 23 . 24 ! 25 ' 26 ! 27 ' 28 ! 29 Much of this work fo- 
cuses on modelling the conditions under which condensa- 
tion can occur, both in equilibrium , 17 i 20 i 21 i 22 anc | consid- 
ering the effects of pumping and decay . 23 ' 26 ' 28 ' 29 Possible 
signatures include the nonlinear relation of emission at 



zero momentum to pumping power f 



27.28.2!) 



changes to the 



linewidth , 26 ' 27 ' 29 the PL spectrum and the angular dis- 
tribution of radiation ; 21 ' 22 and spontaneous polarisation 
of emitted radiation^ 4 - 

In this paper, we discuss the optical properties of con- 
densed polaritons, focusing our interests on absorption, 
photoluminescence (PL) and resonant Rayleigh scatter- 
ing (RRS). From our study we conclude that RRS, the 
coherent scattering by disorder of polaritons into direc- 



tions other than that of the original probe, represents a 
powerful tool for investigating the condensed phase. We 
will show that signatures of condensation are visible in 
the RRS spectrum, allowing a direct probe of the col- 
lective excitation properties of the polariton condensate. 
In particular, we will show that, above the threshold for 
condensation, strong emission from the collective sound 
mode both above and below the chemical potential can be 
seen in the RRS spectra (see Fig. [I]). In contrast such fea- 
tures are expected to be much harder to observe in usual 
PL emission spectra, where the spectrum is dominated 
by the very strong condensate emission at the chemical 
potential, which is likely to mask these more subtle fea- 
tures. In addition, we will show that the RRS spectra 
directly reflects the disorder averaged density of states of 
excitonic particle-hole like excitations, i.e. bound exci- 
tons coupled to the coherent photon field. 

Resonant Rayleigh scattering depends on disorder to 
scatter polaritons between momentum states. In order to 
carry on our analysis, we introduce a realistic description 
of disorder and analyse its effects on RRS and more gen- 
erally on other optical probes, like absorption and PL. To 
do this, we make use of a quantitatively accurate model 
for exciton disorder ) 30 ' 31 i 32 and numerically evaluate the 
distribution of excitonic energies and oscillator strengths 
associated with a given disorder potential. An accurate 
treatment of exciton disorder for quantum wells in mi- 
crocavities is important because the large ratio of exciton 
mass to photon mass means that those exciton states con- 
tributing most to the thermally populated polaritons are 
strongly localised, i.e. cannot be treated within the first 
Born approximation. Such a treatment shows that, for a 
given exciton energy and momentum, there is still a dis- 
tribution of possible oscillator strengths. For the exciton 
states relevant to polariton formation, the distribution 
of oscillator strengths varies from a narrow Gaussian at 
low energies (below the band edge) to a Porter-Thomas 
distribution at high energies, as is known from previous 
worksi 30 i 31 ' 32 Here, we study the effect of this distribu- 
tion on the many-body physics, and in particular, on 
RRS and PL spectra. 

In this paper, we consider the case relevant to the exist- 
ing CdTe or GaAs microcavities of high-quality quantum 
wells, where the typical excitonic disorder amplitude is 
smaller than the Rabi splitting. In this regime, we can 
make use of the coupled oscillator model or wave vector 
conserving approximation, as explained in Ref. l33l . As 
will be explained in more detail later on, this corresponds 
to approximating the full exciton Green's function by 
its momentum-diagonal part (or cquivalently considering 
its disorder average) , and treating perturbatively the off- 
diagonal terms. In this way, the translational invariance 
broken by the presence of the disorder is restored, polari- 
tons have a defined wave vector and, at high enough den- 
sities, condense in the lowest momentum state. Such a 
treatment implies that excitonic disorder, being on short 
length scales, does not lead to spatial inhomogeneity on 
the length scale associated with the polariton. Extended 
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FIG. 1: (Color online) Contourplot of the disorder aver- 
aged RRS intensity (7 pq (a;)) for |p| = |q| as a function 
of the dimensionless momentum |p|a x and rescaled energy 
2(lj — (J,) /{In, for zero detuning, Rabi splitting Or = 26meV, 
temperature fcsT = 20K, and a disorder strength charac- 
terize by an inverse scattering time 1/r — 1.16meV. (a) 
non-condensed regime (dimensionless density p ~ 0); (b) 
condensed regime (p ~ 7.8 x 10 -3 ); (c) condensed regime 
(p ~ 6.7 x 10~ 2 ). (The parameters chosen for these plots are 
the same as those used later for spectral weight and photo- 
luminescence.) The value of the chemical potential is explic- 
itly marked [horizontal green (gray) line]. While in the non- 
condensed regime, RRS emission is always above the chemical 
potential, in the condensed phase, emission from the collec- 
tive sound mode is seen both above and below the chemical 
potential. 



polaritons are formed from a superposition of many lo- 
calised excitons, therefore recovering translational invari- 
ance at the level of polaritons. At low densities, in the 
non-condensed state, many observable properties can be 
adequately described by the coupled oscillator model and 
so can be found from the mean squared oscillator strength 
at a given energy. In this limit our method recovers the 
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well known results for the inhomogencous broadening 
of the polariton PL t 33 i 34 and — by considering in ad- 
dition the mean fourth power of the oscillator strength 
- the averaged RRS responsej 35 ' 36 i 37 ' 38 However, when 
condensed, there are observable effects associated with 
the full distribution of oscillator strengths. In particular, 
the non-vanishing probability of excitons to have arbi- 
trarily small oscillator strengths has direct consequences 
for optical probes, including both RRS and PL. 

In Ref. [2f| we considered specific aspects of resonant 
Rayleigh scattering arising from the model discussed in 
this paper. Here, we provide and compare further ex- 
perimental probes of condensation, and discuss the un- 
derlying physical mechanisms involved. All these optical 
probes, including RRS, are in addition to a non-resonant 
pumping of the microcavity polaritons. A related prob- 
lem is studied in Ref. [3^. There, the Rayleigh scattering 
of a strong resonant pump is considered (treating disor- 
der perturbatively) , thus the probe is the pump. In this 
paper, in contrast, one has to distinguish two types of 
coherence. The first type is the internal coherence of the 
condensed system, which arises spontaneously following 
non-resonant pumping. The second type is the coherent 
scattering of an external laser probe on disorder, which 
can be strongly modified by the presence of the conden- 
sate. 

The form of the density of states, optical density, and 
distribution of oscillator strengths found from the nu- 
merical calculations are also important in other thermo- 
dynamic properties and probes of polaritons. One such 
example is the calculation of the linear blue-shift vs. den- 
sity in the low density regime. The inclusion of exciton 
states that couple weakly to photons extends the validity 
of the current model to densities beyond those in similar 
treatments ) 17 ! 21 and provides a stronger basis for the use 
of those models in the regimes where they are valid. 

The paper is arranged in the following sections: We 
describe the effect of quantum well disorder on the ex- 
citon oscillator strength in Sec. [TTJ There, we also give 
details of the numerical analysis f Sec. Ill A"|) . which will be 
used in the many-body model introduced in Sec. IIIII to 
describe the thermodynamics of polariton condensation 
fSec. HII A"]) . Section HVl describes the optical probes such 
as spectral weight and photoluminescence (Sec. IIV A[) 
and resonant Rayleigh scattering (Sec. IIV C[) . Conclu- 
sions are collected in Sec. |V] while Appendix [Al explains 
how one can detect RRS using phase sensitive measure- 
ments. 



II. EXCITONS IN DISORDERED QUANTUM 
WELLS 

The problem of an exciton in a disordered quantum 
well has been studied at some length in the last two 
dccadeS i 30 ' 31 i 32 ' 40 i 41 ' 42 i 43 Quantum well disorder can arise 
due to interface and alloy fluctuations and affects the 
properties of the excitonic linewidth and absorption spec- 



trum. Similarly to Refs. I30II31I . we will assume the exter- 
nal disorder potential to be correlated on a length scale 
larger than the exciton Bohr radius. Accordingly, we fac- 
torise the excitonic in-plane relative and centre of mass 
coordinates 

V a (r e ,r h ) ~ <p ls (r)$ a (R) 
<p ls (r) = ySMfe- 2 ^ , 

and assume that the disorder affects only the excitonic 
centre of mass motion $ Q (R), while the internal de- 
grees of freedom can be restricted to the ground state 
hydrogenic state tp ls (r). Here, a x = e/e 2 [i r is the ex- 
citon Bohr radius and [i r the reduced mass (henceforth 
we will set h = 1). Neglecting the transverse degrees 
of freedom related to the confinement of the excitons 
in the QW, the energy associated to the wavefunction 
^> a (r e , Yh) is given by the sum of the relative motion en- 
ergy E x [i.e., the band gap minus the exciton binding 
energy, lZy x = (2fi r a 2 )^ 1 = e 4 fi r /2e 2 ] and the energy 
related to the centre of mass motion e a : 



* Q (R) = e Q $ Q (R) 



(1) 



Here, the effective disorder potential V(R) represents 
the microscopic structural disorder averaged over the 
electron- hole motion^ This can be approximated, e.g., 
with a Gaussian noise correlated on a length scale l c > a x 
with variance equal to a 2 : 

l/lc 

(V(R)V(R r )) = (a 2 i 2 c /L 2 ) ^2 e' iq '( R ~ R ') , 

q 

where L 2 is the quantisation area. It will be convenient, 
later on, to introduce the scattering time 

1 



2nvo 2 l 2 ' 



where v = m x /27r is the two-dimensional density of states 
in the clean limit. 

In two dimensional non-interacting systems, all states 
are localised by the disorder potential. However, the lo- 
calisation length and the character of the excitonic wave- 
function change significantly from below to above the 
band edge E x . Well below the band edge, low-energy 
Lifshitz tail states^ are well localised in deep poten- 
tial minima, with a nodeless (roughly Gaussian) shape. 
These states are rare, because they occur due to large, 
rare, fluctuations of the disorder potential. In contrast, 
higher energy states above the band edge have a fractal- 
like shape with many nodes and can be approximated 
by a random superposition of plane waves with the same 
momentum |p| ~ ^/2m- K (e a — E x ). Here, the localisation 
mechanism is closely related to quantum mechanical in- 
terference effects. 

Accordingly, the change of the shape of the centre 
of mass wavefunctions across the band edge is accom- 
panied by a dramatic change in the excitonic oscillator 
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strength^ The oscillator strength describing coupling of 
a quantum well exciton to light, g a . p , is given by the 
probability amplitude of finding an electron and a hole 
at the same position and with centre of mass momentum 
equal to the photon momentum p and is therefore pro- 
portional to the Fourier transform of the centre of mass 



wave function. $ r 



9a,p 



2nu) 



cL 



^i s (0)$ Q , p , 



(2) 



where, d a b is the dipole matrix element. The dispersion 
for photons in a microcavity of width L w is given by 
w p = i/cJq + (cp) 2 /e, where uiq = 2ttc/ L w yfe, and can 
be approximated, for small momenta, by a parabolic dis- 
persion uip ~ w + p /2m p h, with the photon mass given 
by m ph = 2ny/e/cL w . 

Because there may be many different exciton wave- 
functions corresponding to similar exciton energies, the 
oscillator strength g ap is a random quantity, which varies 
both in phase and magnitude. Considering many disor- 
der realisations, we find a distribution of squared oscilla- 
tor strengths, |g QjP | 2 as a function of the energy e a and 
momentum p. This distribution reflects the statistical 
properties of the centre of mass excitonic wave functions. 

To compare to experimentally relevant observables, we 
introduce the density of states (DoS), 



DoS(e) 



V- 



£a)> 



and the mean squared oscillator strength: 

9 2 {e, IPl) = ^^(E I^pI 2 <^ ~ ^)) , (3) 

where (...) is the average over different disorder real- 
isations. These quantities are related to the excitonic 
optical density by the relation: 



D{e) = DoS(e) ff 2 (e,0) 



(4) 



Experimentally, the excitonic optical density can be mea- 
sured by dividing the PL emission by the excitonic occu- 
pation^ 

Focusing for the moment on average properties, rather 
than on the entire distribution, a simple expression ex- 
ists for the mean squared oscillator strength in the high 
energy limit. In this limit, where the DoS is flat and en- 
ergy independent, DoS(e) ~ u, making use of the Born 
approximation, the squared average oscillator strength 
g 2 (e, |p|) does not depend separately on the disorder po- 
tential correlation length £ c and variance er, but instead 
only on the scattering time r (see Fig. [5]): 



9 2 (e,\p\) 



1 



1 



lx ri2 ( £ _ £p )2 + (1/2r) 2 . 



(5) 



where e p = E x + p 2 /2m x is the free particle dispersion. 
From this form, and the effectively constant density of 



states, one can see that both the mean squared oscillator 
strength and the excitonic optical density are symmet- 
rical. The comparison of this approximate form to the 
numerical simulation is shown in Fig. [21 

In contrast, for energies much below the band edge, 
the specific asymptotic expression of g 2 (e, |p|) depends 
whether the correlation length £ c is smaller or larger than 
the localisation length r$(e). In the white noise limit^ 
£ c <C r$(e) ~ (2m x |e — £y) _1//2 j one can show that the 
centre of mass wavefunction <i> Q (R) can be approximated 
by a Gaussian centred at a randomly distributed site R a , 

$ Q (R)=r $ 1 e-( R - R °) 2 /^ 



ip-R Q -(r 4 ,p) 2 /4 



and therefore giving a squared oscillator strength propor- 
tional to: 



|*a, 



2m x |e - -Ex I 



3 -|p| 2 /(4mx|e-Bx|) 



Thus, here the distribution of squared oscillator strengths 
is very narrow, with a mean square value given by the 
above form. 

Similarly, much theoretical (and numerical) work has 
been done to establish the energy dependence of the den- 
sity of states. In the low energy tail, as before, the spe- 
cific asymptotic form of the DoS depends on the value of 
the correlation length £ c . In the white noise limit, one 
can show that DoS(e) oc |e - E yL fl 2 er XXA \ e - B *\l° 2m * i l, 
while, in the opposite (classical) limit, £ c 3> v'^(s) ~ 
(2£ c / y/2m x \e — E^Q 1 / 2 , instead one has DoS(e) cx (e — 

E K ) 2 e-^- E - )2/2a2 (see, e.g., Ref. S3). In general, for 
the finite values of the disorder correlation length cor- 
responding to typical experiments, the regions in energy 
where one of these two analytical regimes apply are very 
restricted, and therefore a numerical analysis is required. 

Numerical analysis is also essential in order to account 
for the distribution of squared oscillator strengths |g Q , P | 2 
near the band edge. This distribution changes substan- 
tially from low energy states to high energy ones. As we 
discuss more in detail in the next section, for high energy 
states the oscillator strength distribution is governed by 
a Porter-Thomas law, while for Lifshitz tail states the 
distribution follows a narrow Gaussian-like distribution. 
Neither of these distributions apply for energies around 
the band edge. However, it is states near the band edge 
that have the largest optical density, and so the distribu- 
tion of oscillator strengths in this region has a significant 
impact on derived quantities. It thus becomes essential 
to use numerical analysis to find the entire distribution 
of oscillator strengths. 



A. Numerical Analysis 

To solve Eq. ((T|), we exactly diagonalise this problem 
within a finite basis set, using conjugate-gradient min- 
imisation with pre-conditioning of the steepest descent 
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FIG. 2: (Color online) Plot showing the energy dependence 
of the excitonic squared coupling strength |<? Q ,p| 2 to photons 
of momentum |p| = (a) (from Ref. I25I . for comparison) 
and |p| = 6.3 x 10 cm - (b), where all exciton states found 
numerically. Results are taken from 160 different realizations 
of disorder potential, and for (b), coupling strengths from 
eight different photon momenta (p x ,p y ) with the same value 
of |p| are combined. The mean squared averaged oscillator 
strength g 2 (e, |p|) for the same value of momentum [lower 
red (gray) points] and the density of states DoS(e) [upper 
green (gray) points] are also explicitly plotted. Inset: Fit of 
g 2 (e, |p|) to the expression ([5]) with a renormalized energy e p . 



wavevector (for a detailed discussion of this method, see 
Ref. [45|). We find the exact eigenvalues and eigenfunc- 
tions in a finite basis on a grid of 120 x 120 points for 
a system of size L — 1/im, a = 2mcV, £ c = 166A and 
excitonic mass m K /mo = 0.08. These parameters give an 
inverse scattering time of 1/t = 1.16meV. For this choice 
of the grid, one can show that convergence is reached. 
From the evaluated eigenvalues e a and eigenstates <& a , P 
over 160 realisations of the disorder potential, we can de- 
rive the excitonic density of states, the oscillator coupling 
strength and its squared average, Eq. ([3]), which wc plot 
in Figs.O and [31 while the corresponding optical density 
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FIG. 3: (Color online) Contourplot of the mean (over 160 real- 
izations of the disorder potential) squared oscillator strength 
g 2 (e, |p|) versus energy and momentum [or equivalently angle 
9 — tan -1 (c|p|/uo)]. Note that the scale in angle is not linear. 
The free particle dispersion i? x + |p| 2 /2m x [solid green (gray) 
line] , and the a trace of the mean squared oscillator strength 
for two representative values of momenta, |p| = (9 — 0°) 
and |p| = 6.3 x 10 s cm" 1 (9 = 82°) [red (gray) plus symbols] 
are explicitly plotted (cf Fig. [2}. The figure is adapted from 
Ref. m. 



is plotted in Fig. |U 

The lower panel of Fig. shows the squared coupling 
strength |g QjP | 2 versus energy for a fixed value of momen- 
tum, |p| = 6.3 x 10 5 cm _1 , corresponding, for a cavity of 
cj = 1.68cV, to an angle of 6 = tan" 1 (c|p|/wo) = 82°. 
Note that, because of the presence of disorder in the 
quantum well, one photon with a given momentum cou- 
ples with many exciton states with different energies. 
We will see later on therefore that a polariton with a 
given momentum is formed by the superposition of one 
photon state |p) and many exciton states |$ Q ). These 
states are more or less strongly coupled depending on 
the distribution of oscillator strength for that given mo- 
mentum. As Fig. [2] shows, by probing a quantum well 
at a large angle, i.e. with high momenta photons, the 
excited excitons with larger oscillator strength are the 
ones which are almost-delocalised in nature and with a 
many-node fractal-like shape. By taking the average over 
many (160) disorder realisations, the squared average os- 
cillator strength <? 2 (e, |p|) shown in the inset in Fig. [2] 
is well described by the Lorentzian shape predicted by 
the Born approximation of Eq. l[5]). with a fitted width 
of 1/Vfit — 1.2mcV in good agreement with the theo- 
retical value 1/t = 1.16meV. However, the peak of the 
Lorentzian does not coincide with the energy of the clean 
limit, e p , but is renormalised down in value, as can be 
shown by employing a self-consistent Born approxima- 
tion. 

In contrast, as shown in the upper panel of Fig. [21 
for photons with zero momentum, the maximum value of 
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the oscillator strength characterises cxcitonic states be- 
low the band edge which are more localised in nature. 
This can be easily understood by the following qualita- 
tive argument: At very low energies, the excitonic state 
is strongly localised in a deep potential minimum and 
has no nodes. Increasing the energy, at first the locali- 
sation length increases [e.g., r$(e) ~ (2m x |e — S x |) -1 ' 2 
in the white noise limit] and thus increases the oscillator 
strength. However, eventually, the wavefunction starts 
developing nodes and consequently the squared average 
oscillator strength decreases. When p = 0, only the 
high energy side of g 2 (e, 0) can be described by the Born 
limit ([5]) . For the chosen values of the disorder potential 
correlation length £ c and variance <r, an analytical ex- 
pression for the low energy (Lifshitz) tail is not known. 

The crossover from localised to more plane-wave-like 
excitonic states obtained by increasing the value of the 
photon momentum is plotted in Fig. [3l At high mo- 
menta, the maximum of the squared average oscillator 
strength g 2 (e, |p|) shown in Fig. [3] follows the free par- 
ticle dispersion, plus a renormalisation down in energy, 
which decreases for higher photon momenta, i.e. these 
states are described well by free particles, including dis- 
order in the first Born approximation. At low momenta, 
the states which have the stronger oscillator strength are 
effectively localised. The crossover, as seen in the con- 
tour plot, happens at relatively large angle, 8 ~ 55°, be- 
cause of the large ratio of exciton to photon mass: The 
crossover momentum is set by exciton mass, but its con- 
version to an angle depends on effective photon mass. For 
a connected reason, the thermally populated polariton 
states are formed out of strongly localised (i.e. beyond 
first Born approximation) excitonic states: For values of 
temperature and photon mass relevant for experiments, 
thermal population of polaritons extends up to around 
10° in momentum, which as seen in Fig. [3] corresponds 
to exciton states not accurately described by the Born 
approximation. For this reason, in the following we will 
concentrate on the oscillator strength corresponding to 
P^O, | 5q: o| 2 . 

Finally we plot in Fig. Q] the optical density. For the 
chosen values of £ c and a, the optical density shows a 
maximum around e* — E x ~ — 0.94meV below the band 
edge, a full width at half maximum (FWHM) of approx- 
imatively a* ~ 0.94mcV and a clear asymmetry of the 
line-shape. It is well know n 32 ' 43 that, by indicating with 
E c = l/2m x £ 2 the confinement energy of the lowest state 
in a typical potential minimum, the excitonic line shape 
is determined by the ratio E c /a and, for a finite value 
of this, the optical density develops an asymmetry to- 
wards higher energies. For our choices of parameters, 
E c /cr ~ 0.85. Asymmetry of the optical density of quan- 
tum well excitons has been also measured experimentally 
(see, e.g., Ref. H3), by dividing the measured PL by the 
excitonic occupation. 

As mentioned earlier, the full distribution of oscillator 
strengths, and not just its mean squared value will be im- 
portant. It is useful to discuss here some technical details 




e - E x (meV) 

FIG. 4: Optical density Eq. Q versus energy. The maximum 
value is around e* — ~ -0.94meV and the FWHM a* ~ 
0.94meV. 

of how this is extracted. The numerical analysis provides 
the excitonic eigenvalues and eigenstates only within a 
finite interval in energy. (The lowest energy states arise 
from rare potential fluctuations, which would require a 
larger region of space to be sampled; the highest energy 
states have spatial variation on lcngthscalcs finer than 
our grid.) Within the interval of energies found, aver- 
ages are performed making use of the raw data coming 
from 160 realisations of the disorder potential; outside 
this interval, averages are taken by extrapolating the nu- 
merics. In particular, in the low energy Lifshitz tail, as 
the oscillator strengths have a narrow Gaussian distribu- 
tion around its squared averaged value, we approximate 
the distribution of |g a ,o| 2 with a delta function at its ex- 
trapolated value g 2 (e, 0). In the very high energy region, 
instead, we make use of the Porter- Thomas distribution, 

tj i I ,2n cxp[-a;/(2x)] 

V{x = \g afl \ ) = -== — , 

V in xx 

where x = g 2 (e,0), again extrapolating the fitted value 
for the squared averaged oscillator strength. In addition, 
we fix the overall scale of |g a ,o| 2 to match integrated opti- 
cal intensity with half the experimentally measured Rabi 
splitting Q,r squared: 

J deD(e) = <J (6) 

This normalisation accounts for factors other than the 
wavefunction in Eq. ([2]). 

Having described the effect of quantum well disorder 
on the energies of excitons and on their dipole coupling 
to the cavity photons, we will next introduce the many- 
body Hamiltonian that will allow us to describe a system 
of thermalized interacting polaritons which, for a given 
density and temperature, can condense to a supcrfluid 
phase. 
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III. THE POLARITON MODEL 

The Hamiltonian describing a model of excitons with 
random energies e a , corresponding to the centre of mass 
cigcnstatc from Eq. |T]), dipole coupled to a cavity mode 
can be written as 

h = y ( b ^ ba + aaa *) + w p^p^p 

a p 

+ ~m Y (g a ,pi/j P bia a + h.c.) . (7) 

* a p 

The description of exciton-exciton interaction in this 
model is somewhat simplified: it imposes the constraint 
that two excitons cannot occupy the same energy level 
E a , because of their composite fermionic nature. This 
constraint is implemented by representing each energy 
level [where energy levels are distributed according to 
the density of states, DoS(e)] by two possible states, re- 
spectively |g.s.) and |ex.). These two states are then rep- 
resented by occupation of two fermionic states, so that 
the ground state |g.s.) = a|jO), and an excitonic state 
|cx.) = b' a \0) = 6J l ,a Q( |g.s.}. Imposing a constraint on 
total fermion occupancy, b^ba + o]jfl a = 1, eliminates 
the unphysical states |0) and at b|jO). In this way sat- 
uration effects are introduced to all orders. Note that 
these fermionic operators aj, cannot be written as a lin- 
ear combination of free electron and hole operators - this 
is a result of having already included the intra-exciton 
Coulomb interaction. 

The transition matrix element between a photon state 
-0 P |g.s.) and an exciton state, (ex.| = (g.s.la^&os, is given 
by <?a, P from Eq. ([2]) and evaluated numerically, as ex- 
plained in the previous section. It is convenient to rescale 
this coupling according to 

\/TZy x m x /2ir , 

where N = lZy x L 2 m x /2TT is the inverse level spacing 
L 2 m x /2ir, measured in units of the excitonic Rydberg 
energy. This corresponds to measuring the density of 
particles per Bohr radius squared. Using these units, we 
may write the total number of excitations 

a p 

n a dimcnsionless form by introducing a dimcnsionlcss 
density of particles p = (N)/N, or cquivalcntly p = 
((N)/L 2 )a 2 .4:7rp/m x , where (N)/L 2 is the physical areal 
density of particles. 

In the limit of vanishing density, our occupancy con- 
straint has no effect, and so this description is equivalent 
to one treating excitons as bosons. As will be shown 
later, in Sec. HVBl for low densities, our occupancy con- 
straint is also equivalent to a model of excitons as bosons, 
but with some effective saturation interaction^ Thus, as 
in the bosonic case, saturation effects are linea r 46 ' 47 and 




e-E x (meV) 



FIG. 5: (Color online) Normalized density of states DoS(e)/v 
and thermal occupation factor [red (gray) line starting at top 
left] for E x - n = 5meV and k B T = 20K. 



the blue-shift increases linearly in density. There is how- 
ever an important difference to those previous works. In 
contrast to the clean systems studied there, our effective 
saturation interaction has a coefficient related to disor- 
der and temperature rather than to the Rabi splitting 
and Bohr radius only. However, increasing the density, 
saturation effects, and thus the underlying fermionic de- 
grees of freedom become more and more important, and 
the conventional linear treatment becomes inadequate. 
In our treatment saturation effects are introduced ex- 
actly and taken into account to all orders. Our model 
is accurate as far as the occupancy is restricted to the 
strongly localised Lifshitz states in the tail below the 
band edge (see Fig. [S]). All the results described here 
respect this limit. One can moreover show that it is this 
density regime which is relevant for on-going experiments 
in CdTepii Higher energy states, beyond the Lifshitz tail, 
contribute to the optical response, however they are only 
weakly occupied, and so the above limit remains valid. 

What the model Eq. Q does not include is the 
Coulomb interaction between excitons with different en- 
ergies e Q , i.e., at low densities, on different localisation 
sites R Q . This contribution is expected to be small in 
the low density regime. This low density regime will be 
observed all the way through this paper. At the same 
time we are not including double-occupancy of a single 
exciton energy level, which could be important at higher 
densities. 

In the following, we will make extensive use of the nu- 
merical analysis of oscillator strength distributions ex- 
plained in the previous section. Let us note here again 
that in the non-condensed regime, the full distribution 
is unnecessary - the excitonic optical density is sufficient 
to derive the polariton dispersion and its inhomogeneous 
broadening. However, in the condensed phase, and in 
order to correctly describe resonant Raylcigh scattering 
mechanisms, we will see that it is of fundamental im- 
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portance to consider the entire distribution of oscillator 
strength. 

We consider a thermal equilibrium system of po- 
laritons, as describing the situation with non-resonant 
pumping and strong thcrmalization processes. This limit 
has been recently demonstrated to be accessible in on- 
going experiments! 13 ' 14 The total density of polaritons 
is fixed by introducing a chemical potential fi. Mak- 
ing use of standard path integral techniques and of the 
grand canonical ensemble, H — liN, we integrate over 
the fermionic fields, and thus write the partition func- 
tion Z = tr[e-^-^], with (3 = l/k B T, in terms of 
the following imaginary time action: 



S[i/>] 



/ dr V i/>* (d T + dp) V'p - Tr In G~ 
Jo 



(9) 



Here, G 1 is the, energy level diagonal, inverse single- 
particle Green's function: 



(10) 



r -i = 9 T +e a /2 Ep5a,P^ P /V^ 
a \E p gt, p r p /VN d T -e a /2 



where s a = e a — /i and w F 



/i are respectively the 



excitonic and photonic energies measured with respect to 
the chemical potential. Within this path integral formu- 
lation, it is possible to show 4 ^ that the constraint on the 
site occupation can be taken into account elegantly by 
shifting the fermionic Matsubara frequencies according 
to 

e„ = (27i + 1)tt//3 h-> e„ = (2n + 3/2)tt//3 , 
which we will assume from here on in. 



A. Mean-Field Phase Diagram 

We will see later in section IIVI how the formation of 
polaritons, the result of strong coupling between photons 
and excitons, influences the optical response and in par- 
ticular how the optical response is affected when the sys- 
tem enters a condensed phase. This section will discuss 
the mean-field equations, and resultant phase diagram; 
however first we wish to stress an important point. As 
observed previously, the relevant polariton states form 
from the mixing of a photon with a given momentum 
to many localised excitonic states with different energies. 
The strength of this mixing is determined by the dis- 
tribution of oscillator strengths. Although the quantum 
well excitons may be disorder-localised, for a weak disor- 
der potential, with a < flu, and in the absence of strong 
photon disorder (which in some cases can also be rel- 
evan t 38 ' 49 ) the resultant polaritons are delocalised. As 
such, polaritons will be described by a momentum quan- 
tum number, and condense, as in the usual picture, in 
the lowest momentum state. 

The static and uniform minimum V-'p( r ) — V"^p,o °f the 
action @ can be found from the equation of motion, or 



saddle-point equation, which has to be solved together 
with the mean-field equation for the total number of ex- 
citations (©: 



UJ = 



1 ,2 1 V- 

—W H > 

N N ^ 



| 2 tanh/3£' a 
2E a 

e a tanh f3E a 



4£ Q 



(11) 

(12) 



Here, the spectrum of particle-hole excitations, i.e. the 
quasi-particle spectrum discussed later, is controlled by 
the energy of an exciton in presence of a coherent photon 
field ip, and is given by: 



E a = yJ(i a /2)* + \g a ,o\ 2 4 2 /N . (13) 

In the following, we will solve simultaneously these 
coupled equations, both above (tp = 0) and below (ip ^ 0) 
the critical temperature for condensation. It is interest- 
ing to compare these with similar equations which have 
been already considered, e.g., in the BCS theory of super- 
conductivity^ or in models describing gases of fermionic 
atoms close to a Feshbach resonance^ 1 - Considering par- 
ticularly the case of Feshbach resonance, the density of 
particles is typically kept fixed, while the interaction 
strength can be varied externally to span across a BEC 
condensate of tightly bound fermionic molecules and a 
BCS state of loosely bound fermionic pairs. In our case, 
two main important differences arise. Firstly, density is 
the parameter which is changed. Secondly, the coupling 
strength of excitons with cavity light depends on the ex- 
citonic energy. There is a wide variety of behaviour that 
can come from these coupled equations pip and (|12|1 . 
controlled by the distribution of oscillator strengths and 
by changing the density of states. 

By solving the 'gap equation' Eq. (jTTJ) in the limit 
tp — > 0, let us concentrate for the moment on the critical 
temperature T c as a function of the chemical potential 
fx, forgetting about the corresponding excitation density. 
At the phase boundary, where ip = 0, in contrast to in- 
side the condensed phase, the squared oscillator strength 
appears only linearly in Eq. pip , this expression there- 
fore involving only the optical density D(e). It is clear 
that the standard BCS result can be thought of as a 
limiting case where the effective optical density is much 
wider than the temperature^ In the opposite limit, of a 
narrow delta-like optical density, one can solve the gap 
equation analytically^ 

By taking into account the density equation (fT2"]) . the 
critical temperature can be expressed as a function of 
the excitation density rather than as a function of the 
chemical potential. In the standard BCS theory there is 
no modulation of the interaction strength and therefore 
no difference between optical density D(e) and density 
of states DoS(e). In our system instead one has to intro- 
duce separately both quantities, and further a realistic 
description of the DoS is required in order to measure the 
density in physical units. In previous treatments; 1 ^ where 
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FIG. 6: (Color online) Phase diagram for the dimension- 
less critical temperature 2ksT c /Or versus the dimensionless 
density p for effective zero detuning S — u>o — £* = OmeV 
(e* is the energy at which the excitonic optical density has 
its maximum, and so as seen in Fig. [3J S = implies 
wo — -E x = — 0.94meV) (solid black) and for positive detuning 
S = +6meV (wo — E x = +5.06meV) (dotted-dashed black) 
and Or = 26meV. The mean-field boundaries for the two dif- 
ferent detunings are cut-off by the expected linear dependence 
of the critical temperature, as indicated at extremely low den- 
sities [orange (gray) solid and dotted-dashed lines]. The hor- 
izontal dashed line marks the temperature of ksT = 20K 
which will be used in later figures. 



critical temperature varies linearly with density, and is 
inversely proportional to the polariton mass: By increas- 
ing the detuning, the polariton mass also increases. 

Having discussed the thermodynamics via the mean- 
field equations, we now turn in the rest of this paper to 
the analysis of the optical responses and probes of micro- 
cavities, both below and above threshold for condensa- 
tion, with particular emphasis on the resonant Raylcigh 
scattering emission. Such optical responses are described 
by considering fluctuations about the mean-field theory 
outlined in this section. The spectrum of these fluctu- 
ations, and the various optical response functions that 
these fluctuations describe is the subject of Sec. II VI 

As well as describing the optical response, fluctuations 
about the mean field theory can also be important in de- 
scribing, in the extremely low density limit, corrections 
to the mean field thermodynamics. Since the subject of 
this paper is optical probes, we will not discuss the sub- 
ject of fluctuation corrections to the critical temperature 
here, as they have been discussed elsewhere^. 

For low densities, a bosonic form can be recovered from 
our Hamiltonian by using the Holstcin-Primakoff trans- 
formation to express the two-level systems in terms of 
bosonic operators, and so this low density limit can of 
course be treated within our model. A smooth crossover 
between the low density (fluctuation dominated), and 
the higher density (long-range interaction dominated) 
regimes can be found by considering fluctuation correc- 
tions to the density, as was studied in Ref l2ll. 



the DoS has been attributed an arbitrary Gaussian-like 
shape (and moreover where the coupling strength was 
kept fixed), the density was measured in arbitrary units. 
Note however that, the fact that the density of states 
asymptotes to a constant at large densities, and the pres- 
ence of those exciton states that couple only weakly to 
light, is responsible only for the asymptotic form of the 
phase boundary at large densities (well beyond the range 
shown in Fig. [5] and beyond the validity of our model). 
At smaller densities, such states have a little effect on 
the phase boundary. In order to calculate the phase 
boundary using our numerically calculated distributions, 
we will solve Eq. (fTTj) and Eq. (fT2)) numerically. How 
the disorder averaged sums in these equations are to be 
evaluated has been discussed Sec. Ill Al 

The critical temperature T c at which the system, for a 
given density, goes from a non-condensed phase to a con- 
densed phase, is plotted in Fig. [5] for two different values 
of the detuning 5 = ujq — e*, where e* is the energy at 
which the excitonic optical density has its maximum (see 
Fig. 2}. It is clear that increasing the detuning of the 
cavity mode above the excitonic one decreases the over- 
all photonic fraction and therefore decreases the value 
of the critical temperature at a given density of excita- 
tions. Conversely, if one fixes temperature, the critical 
density for condensation is higher for larger detuning, 
as Fig. [5] clearly shows. This dependence is easy to un- 
derstand in the extremely low density limit, where the 



IV. OPTICAL PROBES 

The optical absorption and emission spectra of a mi- 
crocavity can be derived by first finding the Green's func- 
tion describing photon propagation. The non-condensate 
part of these response functions (i.e. away from zero mo- 
mentum) may be found in practice by considering the 
Green's function for fluctuations about the mean-field so- 
lution, V P , Wh = V"W<W,o + Sip P ,u h , where u> h = 2nh//3 
is a bosonic Matsubara frequencies, and expanding the 
action §§§ up to quadratic terms: 




This gives the Matsubara inverse photonic Green's func- 
tion: 

Physical response functions can be found by analytic con- 
tinuation of the imaginary time Matsubara Green's func- 
tion to real tiniest The matrix elements of G~* (loh) can 
be expressed in terms of the bare photon energy uj p and 
the excitonic quasi-particle energy E a as follows, 
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sr(i)/ \ r ,■ 1 V- * tanh^gg) z^g a /2 - £g - (e a /2) 2 1 ^ 

^pqK) = 5 P,q («»>/» + W p ) + — 2^ <7 a , p <7a,q ^ ^2 + 4 ff2 <W,° Jf 1^ 7aj., 



^pq(^) = ^ XI l5a : o| 2 5a :P 5a 



a 

V> 2 tanh(/3£; Q ) 



A^ 



r 



(14) 
(15) 



where 



7a : pq = /5|ff a ,o| 2 5a,p5'a, 



V> 2 sech 2 (/3£ a ) 



A^ 



4£ 2 



It is useful to decompose the photonic Green's function 
into its momentum diagonal and off-diagonal contribu- 
tions: 



pp 



pq 



(16) 



By treating the off-diagonal terms perturbatively, the 
translational invariance can be recovered and polariton 
cigenstates can be labelled by momentum vectors. In the 
next section, we will focus on the diagonal terms, which 
characterise both the spectral weight and the photolu- 
minescence emission. Section IIV CI is then dedicated to 
RRS, for which we will see that the off-diagonal terms 
are necessary in order to describe the scattering of an 
incident photon (via the excitonic component of the mi- 
crocavity polariton mode) into directions other than its 
original direction^ 



A. Spectral Weight and Photoluminescence 

Secondary emission from a semiconductor microcavity 
after optical excitation is the source of both incoherent 
PL and coherent RRS. At short times, this emission is 
dominated by RRS, the coherent scattering from disor- 
der, and so is at the energy of, and coherent with, the 
incident radiation. At longer times, phonon and particle- 
particle scattering destroy coherence, and redistribute 
the energy, leading to a quasi-equilibrium distribution of 
energies, and thus the incoherent PL emission intensity, 



P(w,p) = n B {uj)W(w,p) 



(17) 



is given by the Bose occupation factor ub(uj) times the 
spectral weight: 



W(u,p) 



2 



(18) 



The spectral weight can be interpreted as an absorption 
coefficient* 2 ^ (the probability to absorb a photon minus 
the probability to emit a photon) , where negative values 
of W(u), p) represent gain. In contrast, the PL P(w, p) 
is always positive. 

In calculating the PL, it is convenient to make an 
approximation by neglecting multiple polariton scatter- 
ing, while still including the effects of exciton-disordcr 



scattering to all orders. This is discussed in Ref. [33j, 
where comparison between this approximation and ex- 
act numerical calculations show this approximation to 
be remarkably good. Physically, this is a good limit to 
consider because the typical exciton-disorder scattering 
times are very short compared to the inverse frequen- 
cies considered in PL. As PL depends on the momen- 
tum diagonal part of the photon Green's function, ne- 
glect of multiple scattering means in practice averaging 
over disorder realisations at the level of the inverse pho- 
ton Green's function, Gpqi^h)- Since off-diagonal terms 
in the inverse Green's function, Eq. (Til)]) , break transla- 
tional invariance, they average to zero, and so scatter- 
ing between different photon momentum states can thus 
be neglected. The off-diagonal terms neglected here will 
however play a crucial role in the case of RRS response, 
as discussed below. 

The spectral weight calculated from this formula, along 
with the quasi-particle density of states is shown in Fig. [7] 
We first discuss this in the non-condensed case, where 

(2) 

tp = 0, and so Kpq(u>h) — 0, and the photon's Green 
function becomes diagonal in particlc-holc space, and 
simplifies to: 



deD(e) 



tanh(/3e/2) 



lLO h 



This expression describes the coupling of one harmonic 
oscillator (the photon mode) to many harmonic oscilla- 
tors (the exciton modes). In the limit of small density, 
the chemical potential is far below all exciton modes, and 
so tanh(/3e/2) ~ 1. In this limit, the result is identical to 
a bosonic description of excitons sometimes also called 
the linear dispersion model. The underlying fcrmionic 
structure appears as a reduction of the effective exciton- 
photon coupling, due to saturation effects described by 
the tanh(/3e/2) term, and discussed in more detail below. 

From this linear dispersion model, describing one os- 
cillator coupled to many, there are in general two broad- 
ened modes at high and low energies - here these are the 
lower polariton (LP) and upper polariton (UP) - and a 
continuum of modes associated with the exciton optical 
density. However, when considering the corresponding 
spectral weight, modes are weighted by their photonic 
component, so one primarily sees the LP and UP modes, 
and only weak emission near the bare excitonic states 
- the excitonic 'dark' states. These three features are 
clearly visible in Fig. and have been previously pre- 
dicted by Houdre et al^ making use of a simplified and 
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FIG. 7: (Color online) Left column: Contourplot of the spec- 
tral weight W(u,p) as a function of the dimensionless mo- 
mentum |p|a x and rescaled energy 2 (w — fij/Qn, for zero de- 
tuning 5 = ujo-e* = (wo-Sx = -0.94meV), Qr = 26meV, 
fcsT = 20K and photon mass = 0.01. Right col- 

umn: plot of the quasi-particle DoS for the same choice of pa- 
rameters as the left column, respectively: (a) non-condensed 
(p ~ 0, E^ — p ~ 27meV) (the bare exciton and photon disper- 
sions [blue (dark gray) dotted line] and the upper and lower 
polariton dispersions [red (gray) solid line] obtained from the 
effective coupled oscillator model are shown for comparison); 
(b) condensed (p ~ 7.8 x 10~ 3 , E x — p ~ llmeV); (c) con- 
densed (p ~ 6.7 x 10~ 2 , E x - (j, ~ 7meV). 



exactly solvable model. Note that in the limit of an in- 
finitely narrow spectral width, these dark exciton states 
become entirely dark, having vanishing photon compo- 
nent. The presence of 'dark' exciton states, coexisting 
with strong coupling polaritons can explain the simulta- 
neous observation of large Rabi splittings and long decay 
times seen in some experiments^ 

When the Rabi splitting is substantially larger than 
the exciton inhomogeneous broadening, (in our case we 
have fin = 26meV and the FWHM of the optical density 
a* = 0.94mcV), there is a substantial difference between 
broadening of lower and upper polaritons. This is be- 



cause the high energy tail of the optical density decays as 
a power law, while the low energy Lifshitz tail decays ex- 
ponentially. Thus, the optical density has a larger value 
at the UP mode than at the LP mode, giving a larger 
broadening compared to the almost vanishing width of 
the LP. This description of the polariton linewidth due to 
the excitonic inhomogeneous broadening coincides with 
that of Whittaker— , and it has been well tested experi- 
mentally^ 

In addition, the location of the LP and UP can be 
found by making use of an effective two-oscillator model, 
i.e. assuming a narrow delta-like optical density, 

D{e) ~ (^fj 8{e - s*) , 

where e* is the location of the maximum optical density 
and effectively, the exciton energy. In this case the sys- 
tem reduces to two coupled oscillators, giving unbroad- 
ened LP and UP poles at 

Slrup = ± \ v (*p - e* f + n 2 R , (19) 

2 

where , fl R = Ojjtanh(/3e*/2) is the reduced Rabi split- 
ting, due to saturation effects at higher densities. The 
reduction of fl R splitting thus translates directly into a 
reduction of the LP - UP splitting, and thus a blue-shift 
of the LP, which at small densities can be shown to be 
linear. For comparison, the results of this formula are 
shown by the solid (red) lines in panel (a) of Fig. [JJ In 
Sec. IIVBI we will discuss the calculation of this linear 
blue-shift in the low density regime, by evaluating SIr as 
a function of density, making use of the full density of 
states. 

Let us now turn to the signatures of condensation as 
seen in the spectral weight, and thus in the PL emission 
(but which are most probably masked in the PL emission 
by strong emission from the condensate mode). When 
condensed, the polariton modes are replaced by new col- 
lective modes : 21 ' 25 The lower polariton becomes the lin- 
ear Goldstone mode, and two branches appear below the 
chemical potential. The appearance of new excitation 
branches below the chemical potential is generic to con- 
densation^ however the experiments required to probe 
these modes are not easy in other Bose condensed sys- 
tems, such as atomic gases. For this reason, let us briefly 
discuss the physical origin of these new branches, and the 
reason they may be observed in optical response of po- 
lariton systems. The Bogoliubov spectrum arises because 
of the possibility of processes that spontaneously cither 
create or destroy two non-condensed particles. (Such pro- 
cesses arise due to scattering from or to the condensate.) 
As a result of these processes, there is mixing between 
the propagation of an extra particle, or propagation of a 
missing particle (i.e. a hole). (Such a language of par- 
ticle and hole refers to the normal state quasi-particlcs, 
in the current case the normal state polaritons.) In the 
normal state, one can separately calculate the spectral 
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FIG. 8: (Color online) Contourplot of the spectral weight 
W(u>, 0) for p = as a function of the detuning S = loq— e* and 
energy uj — E x , for a fixed and very low value of density (p ~ 0, 
E x - p, ~ 27meV) and k B T = 20K (2fesT/fi fl = 0.13). The 
bare exciton and photon dispersions [blue (dark gray) dotted 
line] and the upper and lower polariton dispersions [red (gray) 
solid line] obtained from the effective coupled oscillator model 
are explicitly shown. 



weight of particle excitations, which have weight only 
above the chemical potential, and hole excitations which 
have weight only below the chemical potential. When 
condensed, as the "particle" and "hole" spectral weights 
become mixed, this mixing can lead to spectral weight 
below the chemical potential that is associated with par- 
ticle propagation. To observe this weight it is however 
necessary to have the ability to inject a particle which 
is not a quasi-particle of the condensed system, i.e. not 
a Bogoliubov quasiparticle. In atomic experiments this 
is hard to achieve, but for polaritons can be naturally 
achieved by injecting a photon. These new branches be- 
low the chemical potential are seen as optical gain in the 
spectral weight (see Fig. [7]). Note that the presence of 
pumping and decay will modify the linear dispersion of 
the Goldstone mode at low momentum, making it diffu- 
sive »2& and that quantisation by disorder may also have 
some effect^ 

The spectral weight also contains information about 
the excitonic quasi-particle DoS of the system. As dis- 
cussed above for the non-condensed case, this is visible 
via the appearance of 'dark' exciton states [see Figs. [5] 
and G0[a)]. When condensed, there is a coherent field 
that modifies the energies of these excitonic quasi-particle 
states, as given in Eq. ([T5]) . and thus modifies their den- 
sity of states. The density of states given by taking the 
energy of such modes is shown by the right column of 
Fig. [Jj and can be compared to the corresponding faint 
features seen in the photon spectral weight in the left 
column. 

The change to the quasi-particle density of states that 
occurs on condensation requires some explanation. As 
mentioned earlier, the change of this spectrum is simi- 
lar to that seen in the BCS theory of superconductivity, 



as such it is surprising that there is no gap in the den- 
sity of states of Fig. [JJ In fact, there would be a gap if 
the distribution of oscillator strengths were replaced with 
the mean square oscillator strength. In that case, there 
would be a branch cut, and a (E — gtp/\/ r N)~ 1 / 2 singu- 
larity in the density of states. However, with the full 
distribution of oscillator strengths, one finds that there 
is always a non- vanishing probability of arbitrarily small 
oscillator strengths. (Note that if the oscillator strength 
were zero, the corresponding state would not contribute 
to any photon response, but for an arbitrarily small cou- 
pling, it has some contribution.) Thus, the contribution 
of these weakly coupled exciton states in effect smoothes 
out the gap. Thus, this system can be vaguely described 
as 'gapless fermion condensation', analogous to 'gaplcss 
superconductivity', but through a mechanism very differ- 
ent to the standard Abrikosov-Gor'kov mechanism con- 
sidered in superconductors. 18 i 50 

As one moves away from the chemical potential, the 
mixing of 'particle' and 'hole' modes described by the 
Bogoliubov spectrum reduces. As a result, far above 
the chemical potential, the modes are much as the un- 
condensed case, and the far below the chemical poten- 
tial, the new modes disappear. To explain this quanti- 
tatively, it is clearer to discuss the case of weakly inter- 
acting Bose gas^ Writing the Bogoliubov mode energy 
as £p = y/ £ p( £ p + 2/it), where \i = gip 2 is the mean-field 
value, the spectral weight is given by: 

W(u,p) = £p + Jt + *> 6{ u - ft,) - £p + Jt ~ 6( U H P ) , 

It is thus clear that at large momenta, where £ p ~ e p + 
the coefficient of the first delta function (modes above the 
chemical potential) will go to one, and the the coefficient 
of the second (modes below) will be suppressed to zero 
(roughly quadratically in energy). 

In contrast to the power law suppression of the spec- 
tral weight of modes far below the chemical potential, 
the PL signal from modes far above zero is suppressed 
exponentially by the thermal occupation of these modes 
[see Eq. (|17p]. while below the chemical potential, there 
is no such decay. However, this discussion neglects emis- 
sion from the condensate mode, which should included 
at zero momentum, as defined in Eq. (|17[) . To sec a no- 
ticeable change of the normal modes requires a relatively 
large condensate density. Due to instrumental broaden- 
ing, the presence of emission from this large condensate 
density might obscure emission from the normal modes, 
so they may only be weakly visible, as shown in Fig. [5] 
As discussed below, RRS may provide a means of escape 
from this problem. 

B. Lower Polariton Blue-Shift 

Before discussing the RRS response, let us briefly dis- 
cuss a point mentioned above - the calculation of the LP 
blue-shift as a function of density, making use of the full 
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FIG. 9: (Color online) Contourplot of the incoherent PL 
P(u, p) as a function of the dimensionless momentum p]a x 
and rescaled energy 2(w — fi)/Sln, for detuning ljo — -Ex = 
-0.94meV, Or = 26meV, k B T = 20K and m ph f2 R a x /2 = 
0.01: (a) non-condensed (p ~ 0, Z? x — p — 27meV); (b) con- 
densed (p ~ 7.8 x 10 -3 , -Ex — p — HmeV); (c) condensed 
(p ~ 6.7 x 10~ 2 , £ x - jU ~ 7meV). 

density of states, and our model of saturation effects ||7J). 
This can be observed in PL experiments by the change of 
blue-shift as a function of intensity of the non-resonant 
pumping. 

In the non condensed regime, the blue-shift of the LP in 
our model is a consequence of the saturation of the (dis- 
ordered) energy levels. Following the narrow band-width 
limit described in Eq. (fl9|) , this saturation blue-shift can 
be found if one has an expression for the chemical po- 
tential at a given temperature and density. In the low 
density limit, the expression for density p as a function 
of \i [Eq. (11211 with ijj = 0] can be inverted in terms of 
elementary functions. As illustrated by Fig. [5j in the 
low density limit, the chemical potential is far below the 



band edge. Thus, there are two significant contributions 
to density: one from low energies, in the tail of the DoS, 
but at large occupation, and one from high energies, in 
the tail of the occupation, but at large DoS, thus: 

f ^ e ~ m/Wp + r =^e-<-">/*- T . (20) 

J-oo ^yx Jo 

Here W p /2 = 0.32meV is the energy which characterises 
the exponential decay of DoS in the tail of the Lif- 
shitz states - this coefficient is extracted by an expo- 
nential fit to our numerical DoS. At a temperature of 
ksT = 20K = 1.72meV, the dominant contribution to 
the density is given by the second integral and the chem- 
ical potential increases logarithmically with the density: 

>~W? UT - < 21 > 

From this expression, and from the coupled oscillator ex- 
pressions of the LP Eq. P^|) , we can explicitly derive the 
reduced Rabi splitting due to saturation effects, 

and the LP blue-shift, which is thus linear in this low- 
density regime: 

5E LP = E LP {p) - E LP {0) ~ Q R n ]— , (22) 

where we have reintroduced here the density n per unit 
volume. Remaining at small densities, but now consid- 
ering the low temperature limit, W p /2 > ksT, then 
the dominant term in Eq. (f2"0"l) is the first, and so 
Eq. (j22|) should be modified by replacing the thermal 
length (mJtBT) , by a characteristic disorder length 

It is instructive to compare this result with that for a 
clean system; in this case blue-shift of the LP has been 
attributed either due to Coulomb interaction^ or to sat- 
uration effects^ (where the expression given here is valid 
only in the dilute limit): 

^Coui ~ "Kyx nal 
SE° at - il. R nal . 

Because the excitonic Rydberg TZy x can be of the same 
order of magnitude as the Rabi splitting f2#, in a clean 
system the two shifts can be expected to be of the same 
order of magnitude. Considering just saturation effects, 
the difference between clean and dirty systems is that, 
at low temperature, in a clean system blue-shift depends 
on exciton number per square Bohr radius, while in a 
disordered system the relevant length is that character- 
istic of the disorder potential, which is in general larger 
than the Bohr radius. Finally, wc wish to observe that 
there is a distinction between Coulomb and saturation ef- 
fects: Coulomb interactions result in a blue-shift of both 
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FIG. 10: (Color online) Contourplot of the spectral weight 
W(u>,0) for p = as a function of the dimensionless density 
p and the rescaled energy u — u>o, for zero detuning (loq — E x — 
-0.94meV) and k B T = 20K (2k B T/fl R = 0.13). The rescaled 
chemical potential (/i — ujq) [green (light gray) dashed], the 
non-condensed lower and upper polariton modes from the two 
coupled oscillator model [blue (dark gray) dotted], and the 
critical density for condensation p c ~ 3.7 x 10 -5 [vertical 
solid red (dark gray) line] are shown for comparison. 



the LP and UP, while saturation leads to a blue-shift of 
the LP and a red shift of the UP, i.e. a collapse of the 
Rabi splitting. Thus, these effects can be experimentally 
distinguished, and their relative magnitudes determined. 

In the condensed regime, the equivalent of the LP is 
the Goldstone linear mode, which by definition starts at 
the chemical potential. Thus, the observed blue-shift is a 
direct observation of chemical potential vs. density. Fig- 
ure [10] shows the variation of the spectral weight at zero 
momentum as a function of density, from which the en- 
ergy of the zero momentum LP and UP modes can be 
extracted. Also shown, for comparison, is the chemical 
potential vs. density (dashed green line). The locking 
of the LP mode to the chemical potential beyond the 
critical density (red solid line) is clearly visible, and the 
behaviour of the LP mode vs density allows the extrac- 
tion of density from an experimental measurement of LP 
blue-shift. 



C. Resonant Rayleigh Scattering 



As described in Refs. [35lr36l . the RRS intensity de- 
scribes the probability to inject a photon into the cav- 
ity with momentum p, and detect a coherent photon at 
q with p / q. The intensity of such a signal may be 
written as: 



/n 



(UJ 



->n i2i 

/pql ' 



I lUJf l — — U) — lT] 



(23) 



Since this involves the off-diagonal momentum space 
components of the Green's function, it is necessary to 
include the off-diagonal momentum components of the 



inverse Green's function Eq. (|16p . However, as discussed 
in Sec. lIVi we will consider only a single polariton scatter- 
ing and again neglect multiple polariton scattering (but 
include all orders of exciton-disorder scattering). Such 
an approximation is reasonable for the same reasons dis- 
cussed at the start of Sec. IIVI 

In the following, we propose using the RRS signal to 
probe the excitation spectrum in the presence of a po- 
lariton condensate. It is therefore necessary to be able to 
separate the RRS signal from the strong photolumines- 
cence that would arise from the equilibrium state with 
large polariton density. Further, since a condensed po- 
lariton system may have a strong nonlinear response to 
an applied probe (i.e. stimulated scattering if the probe 
significantly affects the population of polariton modes), 
it is necessary to use a weak RRS probe. Fortunately, the 
coherent nature of RRS allows exactly this: one can de- 
tect a weak RRS signal by phase sensitive measurement, 
as is discussed in detail in Appendix [AJ Hence, the limit 
on intensity of the probe is provided by the sensitivity 
of the CCD camera, and not by the incoherent photo- 
luminescence background. This therefore allows a probe 
sufficiently weak that only the linear RRS response is 
seen, and nonlinear effects can be avoided. 

When non-condensed, the single polariton scattering 
expansion of the photon Green's function in the off- 
diagonal terms gives 



W^/i) 



|4pW| 2 



l*£, )o H)| s 



(24) 



The factors \K^(u; h )\- 2 appearing here can be inter- 
preted as a filter, allowing a response only when the out- 
going (or cquivalcntly incoming for q — > p) momentum 
has an energy close to the polariton mode at the given en- 
ergy. This means that |p| ~ |q| and so is responsible for 
the ring-shaped RRS signal observed in experiments. 3 ^ 38 
In contrast, the term | Kp^ (tOh ) \ 2 describes scatter- 
ing between momentum states via polariton-exciton- 
polariton scattering. This term has a large variation from 
disorder realization to disorder realization, and is the 
reason for the speckle seen in RRS experiments. This 
speckle, and the disorder averaged RRS intensity are 
shown for comparison in Fig. 111! As the precise speckle 
pattern depends on the precise disorder realization, the 
most we can reasonably do is to describe the statistical 
properties of this speckle, thus the disorder averaged RRS 
signal shown in Fig. [JJ would in experiment represent the 
envelope of the RRS speckle pattern. 

When condensed, the expression for the RRS intensity 
becomes more complicated, but can still be written in 
the form of a filter and scattering part: 



FpSpqFq , 



(25) 
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(Color online) Excitonic speckle pattern 
u h =-w-ir) for a single disorder realization (gray 
dotted line) and its disorder average (black cross symbols) 
versus energy for |p| = |q| = 36° and a 90° azimuthal an- 
gle, p = (0,6.3 x 10 4 )cm _1 and q = (6.3 x 10 4 ,0)cm~\ For 
comparison, we also plot the excitonic optical density [green 
(light gray) plus symbols]. 



where the filter and scattering terms are respectively: 



F n 



Spq — 



I pp 1 



(2)^21 



K {l)o K {l)* K {l)* 

pq pp qq 



-4- K^°* K<-V 

> J1 qp J1 pp J1 qq 



(26) 



pq 1 pp qq pp qq 



As in the non-condensed case, the filter function F p re- 
stricts the allowed incoming and outgoing momenta to 
those for which the normal modes have the injected en- 
ergy. However, as discussed earlier, in the presence of a 
condensate, the normal modes supported are no longer 
the LP and UP, but instead the new Bogoliubov-like 
quasi-modes. 

The decay of spectral weight associated with modes far 
below the chemical potential, as discussed in Sec. IIV Al 
still applies (this can clearly be seen in Fig.[TJ Fig. [7] and 
Fig. [5]). However, since the energy of the RRS signal is 
controlled by the energy of the incident photon, it is pos- 
sible to study modes far above the chemical potential, 
which would have negligible equilibrium occupation, and 
thus negligible weight in the PL signal. Both for this rea- 
son, and because of the ability to discriminate between 
the coherent RRS signal and inchohercnt photolumines- 
cence, RRS provides a very powerful tool to study the 
interesting properties of condensed polariton systems. 



V. CONCLUSIONS 

In this paper we have analysed and compared differ- 
ent optical responses of a condensate of microcavity po- 



laritons. One of the conclusions of this work is that, 
compared to photoluminescence studies, the energy and 
momentum resolved resonant Rayleigh scattering spec- 
tra provides many advantages in probing a condensed 
phase and in studying the associated coherent excitation 
spectrum. In particular, resonant Rayleigh scattering, 
by collecting the signal at a different angle from that at 
which the system is probed, allows one to look directly 
at characteristic features of the condensed phase excita- 
tion spectrum. Such features, like the Goldstone modes 
both above and below the chemical potential, are weak in 
photoluminescence emission and are likely to be masked 
by the strong emission from the condensate in the lowest 
momentum state. In contrast, Rayleigh scattering, be- 
ing coherent with the probe, can be distinguished from 
condensate emission. 

One can ask why the accurate description of excitonic 
disorder considered in this paper is necessary, when the 
microcavity photon introduces a long length scale which 
averages over disorder. One might also compare this sys- 
tem with other excitonic systems, such as double quan- 
tum wellsj^i where high density of dipole-dipole inter- 
acting excitons quickly screens the disorder. 58 ! 59 The an- 
swer to this question is that the large ratio of exciton 
to photon mass means firstly that exciton density re- 
quired for condensation in microcavities is far less than 
in double quantum wells; and secondly that those ex- 
citon states involved in forming the thermally occupied 
polariton modes are strongly localised exciton states (i.e. 
influenced by disorder beyond Born approximation). One 
of the consequences of an accurate treatment of disorder 
in microcavities is that the blue-shift of the lower polari- 
tons due to saturation effects depends on a length scale 
characteristic of disorder and temperature, which can be 
much larger than the Bohr radius, this length playing 
an equivalent role in the clean system. This effect can 
be important in determining the polariton densities from 
the measured blue-shift in current experiments. 

An important consequence of disorder is that the 
exciton-light coupling strength is characterised by a full 
distribution. In the non-condensed phase, the opti- 
cal density alone determines the photoluminescence re- 
sponse, while for resonant Rayleigh scattering the aver- 
aged fourth power of the oscillator strength is required. 
However, in the condensed phase one has to consider the 
entire distribution for each exciton energy. For ener- 
gies close to and above the band edge, there is always 
a non- vanishing probability of an arbitrarily small oscil- 
lator strength. For this reason, in contrast to examples 
like BCS superconductivity, one can show that the quasi- 
particlc spectrum does not have a hard gap. 

In conclusion, we have considered how an accurate 
treatment of disorder on the single-particle excitonic 
level, when elevated to the many body problem of in- 
teracting microcavity polaritons, leads to a variety of in- 
teresting features in various optical responses and probes 
of the condensed phase. 
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APPENDIX A: PHASE SENSITIVE DETECTION 
OF RESONANT RAYLEIGH SIGNAL 

In this appendix we show how the coherent nature 
of the resonant Rayleigh signal allows a weak coherent 
probe to be detected in the presence of strong incoherent 
photoluminescence. The Rayleigh scattering probe is a 
perturbation, H = Ho + V, which in the rotating wave 
approximation may be written as: 



V = A (^ p e mt 



-nit 



This describes a probe at wavevector p, frequency Cl, of 
strength A$. The RRS signal is the coherent scattering 
of this probe to other wavevectors. To isolate the part of 
the emission that is coherent with the probe, one may use 
a homodyne measurement, interfering the emission with 
part of the probe signal. This corresponds to measuring 
the spectrally resolved emission intensity : 



P(w,q)= / dte iut P(t,q) 



(A2) 



where 



P(t,q)=$>^-^) 

71 

x (n\ +A ie -W^] [^(Oj+Ae*] |n) , 

and where A\ is the strength of the homodyne mixing, 
and ^ is a phase delay introduced between the probe and 
the homodyne signal. The states \n) are the eigenstates 
of the system without the probe, given by Hq \ n) = E n \ri); 
F is the free energy, for normalization. 



The principle of phase sensitive detection is to vary the 
phase delay, and to extract the part of P(co, q) which de- 
pends on this phase delay. Since the background PL does 
not depend on the phase delay, this allows one to sepa- 
rate a small but phase dependent signal from a strong 
but phase independent background. 

The part of the signal that depends on phase is: 



P^(t,q)=A 1 J2' 



,0{F-E n ) 



{n\4>l(t)\n)e i 't' + (n|^ q (0)|n>e-^+ m )] . (A3) 

Standard first order time-dependent perturbation theory, 
using the perturbation in Eq. (|A1[) yields: 



(n\^(t)\n) = A J2 
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(A4) 



(Al) The term written [. . .] is similar to the first term, but 
with O i — > — 17 and ip p i— ► tp^. Inserting this signal into 



Eq. (|A3[) and then into Eq. (|A2[) . it can be seen that such 
a term gives a signal at frequency lu = —Cl, and so can 
be clearly separated and ignored. Inserting the first term 
into Eq. (|A3[) , one can write: 



£e«*-*»><ntft(t)|„) 



n.m 
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o 0(E n -E m ) 
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(A5) 



which is the definition of the retarded Green's function^ 

;ing the san 

[31. one can then 



G£ pq ( n ) 



Q}} n I . _ . . By repeating the same 



analysis for the second term in Eq. 
write: 



P*(t,q) = 2A 1 A \G^ p JCl)\ cos(0 + o )e 



-Hit 



(A6) 



where <f>o is the phase of the retarded Green's function. 
Hence, the phase-dependent part of the luminescence is 
given by the off-diagonal in momentum space part of the 
retarded Green's function, i.e. the RRS signal. 



* Present address: Rudolf Peierls Centre for Theoretical 
Physics, 1 Keble Road, Oxford OX1 3NP, UK 

t Present address: Cavendish Laboratory, University of 
Cambridge, Madingley Road, Cambridge CB3 OHE, UK 

1 L. V. Keldysh and Yu. V. Kopaev, Fiz. Tverd. Tela 
(Leningrad) 6, 2791 (1964) [Sov. Phys. Solid State 6, 2219 
(1965)]. 

2 J. J. Hopfield, Phys. Rev. 112, 1555 (1958). 



3 C. Weisbuch, M. Nishioka, A. Ishikawa, and Y. Arakawa, 
Phys. Rev. Lett. 69, 3314 (1992). 

4 Le Si Dang, D. Heger, R. Andre, F. Boeuf, and R. Romes- 
tain, Phys. Rev. Lett. 81, 3920 (1998). 

5 P. Senellart and J. Bloch, Phys. Rev. Lett. 82, 1233 (1999). 

6 Y. Yamamoto, Nature 405, 629 (2000) 

7 P. G. Savvidis, J. J. Baumberg, R. M. Stevenson, M. S. 
Skolnick, D. M. Whittaker, and J. S. Roberts, Phys. Rev. 



17 



Lett. 84, 1547 (2000) 

8 M. Saba, C. Ciuti, J. Bloch, V. Thierry-Mieg, R. Andre, 
Le Si Dang, S. Kundermann, A. Mura, G. Bongiovanni, J. 
L. Staehli, B. Deveaud, Nature 414, 731 (2001). 

9 H. Deng, G. Weihs, C. Santori, J. Bloch, and Y. Ya- 
mamoto, Science 298, 199 (2002); H. Deng, G. Weihs, D. 
Snoke, J. Bloch, and Y. Yamamoto, PNAS 100, 15318 
(2003). 

G. Weihs, H. Deng, D. Snoke, and Y. Yamamoto, Physica 
Status Solidi A, 201, 625 (2004). 

11 M. Richard, J. Kasprzak, R. Romestain, R. Andre, and Le 
Si Dang, Phys. Rev. Lett. 94, 187401 (2005). 

12 M. Richard, J. Kasprzak, R. Andre, R. Romestain, L. S. 
Dang, G. Malpuech, and A. Kavokin, Phys. Rev. B 72, 
201301 (2005). 

13 H. Deng, D. Press, S. Gotzinger, G. S. Solomon, R. Hey, 
K. H. Ploog, Y. Yamamoto, Phys. Rev. Lett. 97, 146402 
(2006). 

14 J. Kasprzak, M. Richard, S. Kundermann, A. Baas, 
P. Jeanbrun, J. Keeling, F. M. Marchetti, M. H. Szyman- 
ska, R. Andre, J. L. Staehli, V. Savona, P. B. Littlewood, 
B. Devaud, and Le Si Dang, Nature, 443, 409 (2006). 

15 F. Tassone, C. Piermarocchi, V. Savona, A. Quattropani, 
and P. Schwendimann, Phys. Rev. B 56, 7554 (1997); F. 
Tassone and Y. Yamamoto, Phys. Rev. B 59, 10830 (1999). 

16 A. I. Tartakovskii, M. Emam-Ismail, R. M. Stevenson, M. 
S. Skolnick, V. N. Astratov, D. M. Whittaker, J. J. Baum- 
berg, and J. S. Roberts, Phys. Rev. B 62, R2283 (2000). 

17 P. R. Eastham and P. B. Littlewood, Solid Sate Commun. 
116, 357 (2000); Phys. Rev. B 64, 235101 (2001). 

18 M. H. Szymanska and P. B. Littlewood, Solid Sate Com- 
mun. 124, 103 (2002); M. H. Szymanska, P. B. Littlewood, 
and B. D. Simons, Phys. Rev. A 68, 013818 (2003). 

19 A. Kavokin and G. Malpuech, Cavity polaritons, in Thin 
Films and Nanostructures vol. 32 (Elsevier, NY, 2003). 

20 F. P. Laussy, G. Malpuech, A. Kavokin, and P. Bigen- 
wald, Phys. Rev. Lett. 93, 016402 (2004); A. Kavokin, G. 
Malpuech, and F. P. Laussy, Phys. Lett. A 306, 187 (2003). 

21 J. Keeling, P. R. Eastham, M. H. Szymanska, P. B. Little- 
wood, Phys. Rev. Lett. 93, 226403 (2004); Phys. Rev. B 
72, 115320 (2005). 

22 F. M. Marchetti, B. D. Simons, and P. B. Littlewood, Phys. 
Rev. B 70, 155327 (2004); F. M. Marchetti, M. H. Szyman- 
ska, P. R. Eastham, B. D. Simons, and P. B. Littlewood, 
Solid State Commun. 134, 111 (2005). 

23 T. D. Doan, Huy Thien Cao, D. B. Tran Thoai, and H. 
Haug, Phys. Rev. B 72, 085301 (2005). 

24 I. Shely kh, F.P. Laussy, A. V . Kavokin, and G. Malpuech, 
preprint cond-mat/0503402 

25 F. M. Marchetti, J. Keeling, M. H. Szymanska, and P. B. 
Littlewood, Phys. Rev. Lett. 96, 066405 (2006). 

26 M. H. Szymanska, J. Keeling and P. B. Littlewood, 
Phys. Rev. Lett. 96, 230602 (2006). 

27 D. Porras and C. Tejedor, Phys. Rev. B 67, 161310(R) 
(2003). 

28 D. Sarchi and V.Savona, AIP Conf. Proc. 893, 401 (2007). 

29 P. Schwendimann and A. Quattropani, Phys. Rev. B 74, 
045324 (2006). 

30 R. Zimmermann, F. Crosse, and E. Runge, Pure and Appl. 
Chem. 69, 1179 (1997); R. Zimmermann and E. Runge, 
Phys. Stat. Sol. (a) 164, 511 (1997). 

31 E. Runge and R. Zimmermann, in Advances in Solid State 
Physics, vol. 38 (Ed. B. Kramer), Vieweg, Braunschweig 
(1998); E. Runge and R. Zimmermann, Phys. Stat. Sol. 



(b) 221, 269 (2000). 

32 E. Runge, "Excitons in Semiconductors Nanostructures", 
in Solid State Physics vol. 57, Academic Press, Amsterdam 
(2002). 

33 D. M. Whittaker, Phys. Rev. Lett. 80, 4791 (1998). 

34 P. Borri, W. Langbein, U. Woggon, J. R. Jensen, and J. 
M. Hvam Phys. Rev. B 63, 035307 (2001). 

35 D. M. Whittaker, Phys. Rev. B 61, R2433 (2000). 

36 A. V. Shchegrov, J. Bloch, D. Birkedal, and J. Shah, Phys. 
Rev. Lett. 84, 3478 (2000). 

37 T. Freixanet, B. Sermage, J. Bloch, J. Y. Marzin, and R. 
Planel, Phys. Rev. B 60, R8509 (1999). 

38 W. Langbein and J. M. Hvam, Phys. Rev. Lett. 88, 047401 
(2002). 

39 I. Carusotto and C. Ciuti, Phys. Rev. Lett. 93, 166401 
(2004); C. Ciuti and I. Carusotto, Phys. Status Solidi B 
242, 2224 (2005). 

40 S. D. Baranovskii and A. L. Efros, Fiz. Tekh. Poluprovodn. 
12, 2233 (1978) [Sov. Phys. Semicond. 12, 1328 (1979)]. 

41 N. N. Ablyazov, M. E. Raikh, and A. L. Efros, Fiz. Tverd. 
Tela (Leningrad) 25, 353 (1983) [Sov. Phys. Solid State 
25, 199 (1983)]. 

42 M. Litinskaia, G. C. La Rocca, and V. M. Agranovich, 
Phys. Rev. B 64, 165316 (2001). 

43 R. Zimmermann, E. Runge, and V. Savona, "Theory of 
resonant secondary emission: Rayleigh scattering versus 
luminescence", in Quantum coherence, correlation and de- 
coherence in semiconductor nanostructures (Ed. T. Taka- 
gahar), Elsevier Science (2003). 

44 I. M. Lifshitz, Adv. Phys. 13, 483 (1964); J. Zittartz and 
J. S. Langer, Phys. Rev. 148, 741 (1966); B. I. Halperin 
and M. Lax, Phys. Rev. 148, 722 (1966). 

45 M. C. Payne, M. P. Teter, D. .C. Allan, T. A. Arias and 
J. D. Joanopoulos, Rev. Mod. Phys. 62, 1045 (1992). 

46 G. Rochat, C. Ciuti, V. Savona, C. Piermarocchi, A. Quat- 
tropani, and P. Schwendimann, Phys. Rev. B 61, 13856 
(2000). 

47 C. Ciuti P. Schwendimann, B. Deveaud, and A. Quat- 
tropani, Phys. Rev. B 62, R4825 (2000). 

48 V. N. Popov and S. A. Fedotov, Zh. Eksp. Teor. Fiz. 94, 
183 (1988) [Sov. Phys. JETP 67, 535 (1988)]. 

49 M. Gurioli, F. Bogani, L. Cavigli, H. Gibbs, G. Khitrova, 
and D. S. Wiersma, Phys. Rev. Lett. 94, 183901 (2005). 

50 A. A. Abrikosov, L. P. Gor'kov and I. E. Dzyaloshinski 
Methods of Quantum Field Theory in Statistical Physics. 
(Dover, New York) 1975. 

51 M. Holland, S.J.J.M.F. Kokkelmans, M.L. Chiofalo, and R. 
Walser, Phys. Rev. Lett. 87, 120406 (2001); E. Timmer- 
mans, K. Furuya, P.W. Milonni, and A.K. Kerman, Phys. 
Lett. A 285, 228 (2001); Y. Ohashi, A. Griffin, Phys. Rev. 
Lett. 89, 130402 (2002). 

52 P. B. Littlewood, P. R. Eastham, J. M. J. Keeling, F. M. 
Marchetti, B. D. Simons, M. H. Szymanska, J. Phys.: Con- 
dens. Matter 16 S3597 (2004). 

53 R. Houdre, R. P. Stanley, and M. Ilegems, Phys. Rev. A 
53, 2711 (1996). 

54 A. V. Mintsev, L. V. Butov, C. Ell, S. Mosor, G. Khitrova 
and H. M. Gibbs, Pis'ma Zh. Eksp. Teor. Fiz. 76, 742 
(2002) [Sov. Phys. JET Lett. 76, 637 (2002)] 

55 L. P. Pitaevskii and S. Stringari, Bose Einstein Condensa- 
tion (Clarendon Press, Oxford, 2003). 

56 A. V. Kavokin. Private Communication. 

57 L. V. Butov, J. Phys.: Condens. Matter 16 R1577 (2004). 



18 



A. L. Ivanov, Europhys. Lett. 59, 586 (2002). 

R. Zimmermann, Solid State Comm. 134, 43 (2005). 

Formally this can be seen by representing the two levels as 



a spin 1/2, and making a Holstein-Primakoff transforma- 
tion of our model. 



